Spin-orbit excitations of quantum wells 
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Confinement asymmetry effects on the phot oabsorpt ion of a quantum well are discussed by means 
of a sum-rules approach using a Hamiltonian including a Rashba spin-orbt coupling. We show that 
while the strength of the excitation is zero when the spin-orbit coupling is neglected, the inclusion 
of the spin-orbit interaction gives rise to a non zero strength and mean excitation energy in the 
far-infrared region. A simple expression for these quantities up to the second order in the Rashba 
parameter was derived. The effect of two-body Coulomb interaction is then studied by means of a 
Quantum Monte Carlo calculation, showing that electron-electron correlations induce only a small 
deviation from the independent particle model result. 

PACS numbers: 



Quantum wells, i.e. strongly confined electrons at 
a semiconducting heterostructure, are intrinsically very 
complex systems, in which the local arrangement of the 
electronic structure can in principle give rise to a large va- 
riety of interesting subtle phenomena. However, the very 
strong confinement and the weak potential felt by the 
electrons allowed for successfully using simplified Hamil- 
tonians in describing the phenomenology of these sys- 
tems. The most celebrated is the two dimensional elec- 
tron gas (2DEG) in the effective mass and dielectric con- 
stant approximation, and its numerous variants. At the 
semiconductor junction, the effective potential is never 
symmetric with respect to the junction plane. This im- 
plies that electrons feel on average the effect of an elec- 
tric field transverse to the plane. In the frame of effective 
Hamiltonians, this fact is accounted for by the inclusion 
of the so-called Rashba term, which simply describes the 
interaction between the spin of an electron moving in the 
plane, and the self induced magnetic field. For TV elec- 
trons the 2DEG Hamiltonian reads: 
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where e and m* are the effective dielectric constant and 
electron mass in the semiconductor, respectively. The 
Rashba interaction I2JI : 
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where is the momentum of the i — th electron, and 
erf and of are the Pauli matrices acting over the spin of 
particle i, has a strength A proportional to the modulus 
of the transverse electric field. Given the formal analogy, 
the Rashba term is usually reported as a spin-orbit (SO) 
interaction. The study of effective SO effects in confined 



electron nanostructures has been the object of many ex- 
perimental and theoretical investigations in the last few 
years. The most interesting fact is that SO links the spin 
and charge dynamics, hence opening the possibility of 
spin control by means of electrical fields. It affects the 
charge transport lln3 and the far-infrared absorption in 
heterostructures fsPfl2| and can lead to spin inversion (13- 



17| and variations of the electric spin precession in a mag- 
netic field fl8M22l ] . It also induces linear magneto-optical 
Kerr or Faraday effects then allowing the measurement 
of the ma gnet ization of materials by means of laser pulse 
excitation|23(. In this sense, the Rashba 24 SO contri- 
bution has received special attention since it was shown 
that its intensity can be externally tuned via the appli- 
cation of a gate voltagejifl Eij . However, it would be 
very important in view of technical applications to have 
a clean and model independent way to measure the ini- 
tial strength of the Rashba SO interaction in the sample, 
in order to quantitatively control the induced effects in 
the physical processes of interest. 

In this work we show that this measurement can be 
achieved by far infrared spectroscopy in the quantum 
well, in absence of external magnetic fields. If the trans- 
verse electric field, either given by the confinement asym- 
metry or external, is absent, then the light absorption 
is zero. On the contrary, when the SO interaction is 
present, photoabsorption occurs with a strength and at 
an energy in the far infrared domain which are linearly 
proportional to the transverse electric field strength. The 
actual experiment would provide, in addition, a very 
stringent test on the applicability of the Rashba model 
to real systems. 

In the following we will consider a quantum well 
modeled by N electrons laterally confined in the xy 
plane, at density p, and interacting via the Hamilto- 
nian H = H2DEG + VR as hba- We make use of effective 
atomic units in order to simplify the notation, defining 



h = e 2 /e 



= 1. The lenght unit therefore is the 



effective Bohr radius ag = age/m* and the energy unit 
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the effective Hartree H* = H ■ ?n*/(m e e 2 ). 

The single particle part of Hamiltonian ([!} can be 
solved analitically[24j to give the following solutions for 
single particle energies and wave functions: 
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where ip = tan^ 1 -A. Note that the wave function ([3]) is 
an eigenstate of a component of the spin which is per- 
pendicular to the direction of k. In the following we will 
address these states as quasi-up and quasi-down single 
particle states respectively. For a given chemical poten- 
tial /i > 0, the two branches of the excitation energy t k± 
, are filled up to the Fermi momenta k^ and k^_, where 



k F ± = y/2np(l T £ 
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(see for example Ref|2j 

We now turn to analyze the transitions induced in the 
quantum well by the interaction with a linearly polarized 
(along y) electromagnetic wave propagating along the x- 
direction, i.e., in the plane of motion of the electrons. 
The corresponding vector potential is A(i) = AsmOj, 
with 9 = tut — qx. The interaction Hamiltonian J • A/c + 
g* hb s • (V x A), where J = v, reads 
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and the velocity operator v y is defined as v v = —i[y, H] = 
p y + \a x . 

This interaction Hamiltonian yields three excitation 
operators: X^Pii a f an< 4 12i °f- The first one does 
not excite the electrons of the system since it commutes 
with the quantum well Hamiltonian. The second opera- 
tor is sub-leading in A with respect to the third operator, 
and is normally negligible. We are then left with the 
excitations induced in the target by the operator 
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Note that for a different geometry, e.g. the one corre- 
sponding to transitions induced in the system by the in- 
teraction with a left-circular polarized electromagnetic 
wave propagating along the z-direction perpendicular to 
the plane of motion of the electrons, the relevant excita- 
tion operators are 



(7) 



where cr* = a x ± ia y . For time- reversal invariant sys- 
tems, like the ones considered here for which (OjS^lO) = 



0, it is possible to show[28| that the strengths and mean 
excitations energies corresponding to the operators ([?]) 
coincide with the ones induced by the operator ©. 
Photoabsorption experiments attempting to measure the 
strength A of the SO interaction could then be performed 
with electromagnetic waves propagating in a direction ei- 
ther parallel or perpendicular to the quantum well. 

The strength and mean excitation energy of the tran- 
sitions induced in the quantum well by the operator© 
can be studied in a convenient way by means of the mo- 
ments nik of the dynamic form factor S(S z ,uj) (which is 
defined as J2 n I < 0\S z \n > \ 2 5(u - w„)): 
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where uj no = E n — E is the excitation energy. 

By using the completeness relation ^ n |n)(n| = 1 and 
the equation H\n) = E n \n), it is possible to write the 
moments (JSj> as the mean values on the ground state of 
commutators [ ] and anticommutators { } of the excita- 
tion operator S z and of the Hamiltonian H. For time- 
reversal-invariant systems, the following sum rules are 
derived for the operator S z : 
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mi = -(0\[S',[H,S']]\0), 
m = U0\{S z ,S z }\0), 



m_i = -(0\[[X,H],X]\0), 



(9) 



where the operator X is the solution of the following 
equation 



[H, X] = S z 



(10) 



The moments can be used to define different mean 
excitation energies, namely 



J5i,_i = \fm~Jm~i 



and 



Ei t o = m 1 /m , 
which satisfy the relations 
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These energies strictly coincide only when the whole 
strength is exhausted by a single state. However, in pres- 
ence of collective states, they can be used to study the 
spreading of the excitation strength |28[ and to bind the 
energy of the collective modes. 

In the single particle picture and to the lower order 
in A, the momenta mi, mo, and m_i can be computed 
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analytically. The expressions are: 

™i = -±{0\V Raahba \0) = ^ 2 (l-±g), 
mo = (\ — \—\ , (14) 

= S; ■ 

The corresponding mean excitation energies (|12l) and 
(jT2j) become: 

£i,o = =0^A(l-ig), 

(15) 

= =^A(l-ig). 

Some comments are in order here. The total strength too 
and the energies (fl~5"T) vanish for A — > 0. As anticipated 
in the introduction, light absorption in the system only 
occurs in presence of SO interaction and the excitation 
strength and energy are linear in A at the main order. 
This fact characterizes this excitations as a genuine effect 
of the well asimmetry, or of the presence of a transverse 
external magnetic field described by the Rashba SO term. 

At the lower order in A the two energies fTS"]) coincide, 
and the higher order corrections are close each other and 
very small (of order A 2 ). As a consequence the excita- 
tion strength is practically concentrated in a single peak. 
This peak is measurable since it lies in the far infrared 
region. In fact taking, as an example, a GaAs quan- 
tum well for which e = 12.4 and to* = 0.067TO e , yclding 
H* = 11.86TOey and a* Q = 97.93i, one gets for E lfi a 
value ranging from 2 to 6 meV, for densities in the range 
(2 13)10 n cm- 2 and A ~ lO^eVcm ~ O.la.u.. 

We want to stress the fact that the excitation is linear 
in A only in absence of an external magnetic field, which 
would immediately make such dependence quadratic, 
therefore strongly affecting the possibility of accurately 
measure the SO strength. We also stress that, since the 
excitation operator S z commutes with the Coulomb in- 
teraction, the result mi = — \ (0| VR as /j;,o|0) for the en- 
ergy weighted sum rule is exact for the total Hamilto- 
nian ([1]). Deviations from the result (fT4"f for mi are due 
to two-body correlations induced in the ground state by 
the Coulomb interaction. Finally we note that result 
(fl~4| for m_i, valid at all orders in A, is exact for the 
single particle model where the Coulomb interaction is 
neglected. 

In order to evaluate the effects of the two-body 
Coulomb interaction on the strength too and the energy 
i?i.o we have performed a Quantum Monte Carlo (QMC) 
calculation. Our recent extension of the Diffusion Monte 
Carlo (DMC) [2j|, in fact, allows for the projection over 
the ground state of the system corresponding to the full 
Hamiltonian ([T]). DMC is a very accurate method to 
solve the Schrocdinger equation, based on projecting out 
of an initial trial wavefunction ^t{R) the ground state 
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TABLE I: Extended DMC ground state estimates for the first 
and second moments of the structure function and for the 
average excitation energy (in H*) for the excitation operator 
S z as a function of the SO strength A. Beside A, the number 
of "quasi- up" (JV+) and "quasi-down" (JV_) states used in the 
trial wave function are reported. 

of the Hamiltonian by means of an importance sampled 
imaginary time propagation: 

*t(R)4>(R,t) = J dR'G{R;R';r)^^r^ T (R')cP(R,0) 

(16) 

In this case R indicates both the space and spin coordi- 
nates of the electrons. The Green's Function for short 
imaginary time can be approximated by: 

G(R,R',At) OC e-VW^aArg 2i4x e -(V Coui Ar-,Eo) 5 

(17) 

where Vcoui is the Coulomb potential, and Eq is a con- 
stant used to preserve the normalization. The factor in- 
cluding Vftashba contains an explicit dependence on the 
spin coordinates of the electrons. This fact demands 
a substantial modification of the standard algorithm, 
which must now take care of the propagation of the spin 
degrees of freedom together with the space coordinates. 
The expectation giving the quantities of interest can then 
be evaluated as stochastic averages. We should mention 
that, as in all applications of standard DMC to many 
Fermion systems, one must introduce an artificial con- 
straint (in this case a fixed-phase constraint) in order 
to avoid the notorious sign problem. The details of the 
algorithm and further discussions can be found in Ref. 

We have evaluated the momenta mi and mo by explic- 
itly calculating the expectation values of eqs. (|9|) on the 
Monte Carlo ground state at r s = 1, corresponding to a 
density p ~ 3 • 10 11 cm -2 . Calculations were performed 
using the trial wave function described in [29(, using 
the combination of quasi-up and quasi-down single 
particle states giving the lowest ground state energy. 
The obtained results are reported in Table Q] and Figure 
[T]for three values of A expressed in effective atomic units. 

In Fig. 1 we compare the A dependence of the excita- 
tion energy predicted by the independent particle model 
(with and without corrections at order A 2 ), and by the 
DMC calculations. As it can be noticed, the deviations 
due to the presence of the electron-electron correlations 
are rather small. At values of A of order of the estimates 
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FIG. 1: (Color online) Average excitation energy estimates 
are reported in H* units. The green line and the black dots 
correspond respectively to the lowest A order estimates (lin- 
ear) and to the exact values for the independent particle 
model. The red squares correspond to the QMC estimates 
in presence of Coulomb interaction. 



for typical quantum wells they are almost negligible, and 
increase with the strength of the transverse electric field. 
This implies that the analytical results ([14"]) and (fT5j) can 
be considered very good estimates for the strength and 
excitation energies of far infrared absorption of polarized 
light in the quantum well. This confirms the fact that 
a good comparison with experiments would support the 
physical content of the 2DEG approximation even for as 
concerns the description of the effects of the well asym- 
metry. 

In summary, we have calculated strength and mean 
excitation energy of the transitions induced in a two- 
dimensional quantum well by a linearly polarized electro- 
magnetic wave propagating in the plane of motion of the 
electrons by means of sum-rules techniques. In absence 
of a transverse electric field there is no absorption of light 
in the system. When the electric field is turned on, and 
the corresponding strength of the Rashba SO coupling is 
non-zero, the predicted strength and mean excitation en- 
ergy turn out to be proportional to the Rashba coupling 
parameter A. Calculations have been performed analyt- 
ically in the single particle picture which neglects the 
Coulomb interaction and takes into account the SO po- 
tential. The effects of the Coulomb interaction have been 
accurately tested by numerical Diffusion Monte Carlo 
simulations including the recently developed extension 
to the treatment of spin-orbit interactions. 
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Wiglaf HPC facility of the Physics Department of the 
University of Trento and on the CINECA HPC under a 
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